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Abstract 

This paper proposes a sensible definition of a deformation metric between 2-dimensional 
surfaces obtained from each other by an area preserving (incompressible) mapping, and 
an algorithm for obtaining this metric, as well as the optimal deformation. 

1 Introduction 

Recently, there is an increasing interest in the analysis of near-rigid deformation within the 
computer vision community, in particular pattern recognition, image segmentation and face 
recognition ([B2K], [EK], [MS]). 

For example, a variety of objects can be represented as point clouds. These can be 
obtained by sampling of the objects in question by points in some canonical Euclidian space. 
One is often presented with the problem of deciding whether two of these point clouds, and/or 
the corresponding underlying objects or manifolds, represent the same geometric structure 
(object recognition and classification). 

To quantify the difference between two such clouds it is natural to construct smooth 
domains out of the samplings, and look for a mapping from one domain to the other which is 
" as close to an isometry as possible" . The minimal deviation of these mappings from a rigid 
deformation (isometry) can stand as a measure of similarity between the original objects. 

However, this task is very difficult from a computational point of view, since the set of all 
mappings between two domains is very large. On the other hand, it is sensible to assume that 
the density of the sample points refiects the true nature of the object. This implies that the 
volume elements associated with these domains, created out of the samplings, are prescribed. 
In particular we may restrict ourselves to an approximation of an isometry which preserve 
the volume of the two domains. 

The object of this paper is to propose a sensible definition of a deformation metric be- 
tween 2-dimensional surfaces obtained from each other by incompressible (in this case, area 
preserving) mapping, and an algorithm for obtaining this metric, as well as the optimal 
deformation. 

In section [2] the problem is formulated for a pair of flat domains. Section [3] introduces the 
analytic conditions for quasi-rigid deformation, for flat domains. In section [J] the problem 
is extended to a pair of embedded surfaces in M'^. In section [5] we use the results of the 
previous sections to propose an algorithm, namely a flow which converges (formally) to a 
quasi rigid deformation. For the convenience of the reader we defer all technical proofs to 
the final section [6l 
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2 A particular example 

Consider, for example, two flat domains Q, Jli C M^, equipped with the Euclidian metric 
e{dx,dy) = dx^ + dy"^. If Q.i is a rigid deformation of O then there exists an isometry 
$ : (Jl, e) — > (r2i,e). Locally, it means 

= = $^-^>j, = onfi . (2.1) 

In general 11 and f^i are not isometric, and there is no mapping verifying (j2.1|) . However, $ 
is, by definition, an isometry between (11,(7$) ^-iicl (lli,e) where 

g^{dx,dy) := { dx dy ) \D^\^ ( ) = I^^I^^^^ + • ^y)dxdy + |«>yp(iy^ . (2.2) 



1$ P / 



I r,$|2 ._ \^x\ ^x-^y 



How can we quantify the deviation of (r2i,e) from an isometric image of (r2,e)? Certainly, 
it is related to the deviation of 5$ from the Euclidian metric. Recall that 5$ is represented 
by a symmetric 2x2 matrix. So, we consider a real valued function defined on the set of 
symmetric 2x2 real matrices S{2] R). Let ha : S{2] M) R, where a is some real parameter 
(see below), verifying 

K{A) = ha{UAU*) for any U G 0(2;R) and A G 5(2; R) , 
K{A) > K{I) = 

where / is the identity 2x2 matrix, and the equality holds if and only if yl = /. Now set 

HJn,^i)= inf I hJ\D^^)dxdy . (2.3) 

$eZ)i//(f7;Ci) 

By this definition we obtain that (11, e) and (Hi, e) are isometric if and only if Ha{^}, Hi) = 0. 

What is a natural choice of /i^? It must be a function of the eigenvalues Aj, i = 1,2 of 
A, hence it depends on only 2 arguments, say tr{A) and det{A). Note that 

k{A) := tr\A) - Met{A) = (Ai + As)^ - 4A1A2 = (Ai - As)^ > , (2.4) 

and A = I if and only if k{A) = and det{A) = 1. This leads us to the natural choice 

ha{A) = k{A) + a {det{A) - if , 

where a > is a parameter. 

In this paper we restrict ourselves to incompressible deformations. This corresponds to 
the choice a = 00 which implies the constraint det{A) = 1. The adaptation of the quasi- 
rigid deformation metric ()2.3[) to the incompressible case is obtained by the constrained 
optimization 

Hoo(n,ni):= inf [ k(\D^^)dxdy (2.5) 

$GO(f7;ni)7n 
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where 0(0; Jli) is the set of ah area preserving diffeomorphisms ^> : i7 — > fii, that is, 

0{VL, VLi) := E DiffiQ; l^i) ; - ■ ^yf = 1 on . (2.6) 

Can we compare the domains il, 0,i using the definition (j2.5p ? Evidently, Hoq{0,; < oo 
if and only if the set 0(0; Oi) is non empty. By its definition, a necessary condition is 

Area{n) = Area(Oi) . (2.7) 

By a theorem of Moser [M] it tm'ns out that condition (j2.7p is also sufficient, under rather 
general conditions. 

Assuming (|2.7p . the existence of a minimizer of (j2.5p is a much more difficult problem. 
By the definition (j2.4p of A;(-) and ()2.6p we pose the following equivalent open problem. 



Open problem: Suppose \2. 7| ). /s i/iere a minimizer of 

inf 

*eO(n;Ci) 



inf ^($) w/iere i?($) := I tr'^ i\D^^) dxdy ? (2. 



3 Main result for flat domain 



The set (|2.6p is, formally, an infinite dimensional manifold. There is an associated right- 
translation on this manifold by the group of area preserving diffeomorphisms 

0(0) := {S G Diff{^) ; \S:,\^\Sy\^ - {S^ ■ Syf = 1 on 0} . (3.1) 

Indeed, 0(0) is a group under composition, and its action on 0(0;0i) from the right is 
defined by 

$ G 0(0; Oi) ^ $ o 5 G 0(0; Oi) . 
0(0) is a formal Lie group. Its Lie algebra is given by 

o(0) := {v G C'°°(0; M^) V • = on O , u • n = on 90 .} . 
Any smooth flow t G 0(0; Oi) is generated by an orbit t v^^^ G o(0) via 

1$W=^W . (3.2) 



The flow (j3.2p can be presented by the Euler equation 

^^(t) _ . = . (3.3) 

dt ^ ' 

Our object is to deflne v^^^ G o(0) for which 11 is strictly decreasing, and is stationary 

if and only if is a critical point of il in the sense described below. 

At this stage it is convenient to consider complex number notation. Here we represent 
(x, 2/) ~ X + iy G Z. Then 

dz = -{dx—idy) ; dz = -{dx+idy) , dz = dx+idy, dz = dx—idy, dzf\d:z ~ —2idxdy . (3.4) 
Lemma l3. II below is elementary from the divergence theorem and the deflnition (j3.4ll : 
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Lemma 3.1. // O is simply connected then for any v G o(r2) there exists a function ip € 
C°°(r2) so that Ip = on dO, and v = {l/2){—tl'y,ipx) on In particular, ~ f i +iv2 = idjip 
and 

idzipdz = -—\dz\ on 317 
on 

where dtp/dn is the (real) outward normal derivative of ip on dQ. 
For ^> € 0{n; ^i), let the Hopf function [H] 

f^{z,z):=\%\^-\^^\^ + 2i{^,-%) . (3.5) 



The main result is: 

Theorem 1. Let <1>° G 0{0,;0,i). Let <i>W be a flow \3. 3\) where ~ idzijj^^^ satisfying 
e C~(0), V^*^ =0 ondn for any t G M. Then $W G 0(0; Oi) for any t G M and 

^^H (<!>(*)) = 4am ^ f^(t)d^i,dxdy . (3.6) 

Now, we are in a position to define quasi-rigid deformation as a critical point of the 
functional (j2.8p on the constraint manifold (j2.6p . 

Definition 3.1. ^ mapping <I> G 0(f7i; ^2) is quasi-rigid if and only if dH{^)/dt = /or any 
^ g C°°(r2) verifying ip = on dQ. 

Proposition 3.1. ^ : ^ Qi is a quasi rigid deformation if and only if {2.1^ and the 
following conditions hold: 

i) 3m (a|/$) =0 on^. 

ii) 3m (/$g) =0 on dn. . 

Remark 3.1. It is interesting to compare the conditions of Proposition to the case of 
harmonic maps between two Riemannian surfaces. In that case the function /$ is holomor- 
phic, that is, djf^ = 0. This follows from the stationarity of the Dirichlet functional H with 
respect to all parameterizations of the domain il. In the case under consideration, the func- 
tional H is constrained to the set 0(0,), and the criticality condition (i) of the Proposition is 
weaker. Note also that (ii) can also be written as 3m{f^dz'^) = on dQ, while f^dz"^ is the 
Hopf differential (see Definition 1.3.10 of [H]). 

4 Generalizations 

Here we generalize the results of section [3l Instead of flat domains (7, l^i, we consider a pair 
of smooth, compact surfaces S, Si embedded in M^. 

Assume that S, Si are diffeomorphic to a canonical domain A. For simplicity we may 
think about the case A := {x^ + < 1} C M^. We shall later comment about the case where 
0,J7i are embedded manifolds without boundary (e.g. A is the unit sphere S^). 
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Let 

X : A ^ S = X(A) C , Y : A ^ Y(A) = Si C , 
difFeomorphisms from A to S and Si (respectively), then the mapping 

YoX-^iS^Si (4.1) 

describes a diffeomorphism between the surfaces S and Si. 

In order to investigate the deviation of Y o X~^ from an isometry, we use the mappings 
X and Y to pull back the metrics from S and Si respectively, to A. The geometry of S is 
pulled back to A by the parameterization X as 

gy^{dx,dy) := { dx dy ) G'xl(^^y) ^ = + 2(Xa; • Xy)dxdy + W'Xyfdy'^ . 

(4.2) 

where the standard inner product (•, •), H^lP := {v,v) in is used, and 

^ I _( ll^xlP • Xy) \ 

^^l(-^)^-^(x..x,) !|x,f )^^^^^ ■ 

The corresponding area element pulled back from S to A via X is 

|Xx- A Xyl dxdy = det (Gx) dxdy . 

Similar expressions hold for the metric induced from Si on A via Y. 

Given such reference parameterizations X and Y, the set of all diffeomorphisms from S 
to Si can be represented by the set of all diffeomorphisms $ : A — > A as follows: 

Yo$-ioX"^ :S^Si , ^(£Diff{A). (4.3) 

Given $ = {(pi,(p2) G Diff{A) we inquire the deviation of the mapping (|4.3|) from an 
isometry. The pull-back of the metric from S to A by X o $ := X* : A — > S is given by 

Gx* \icc,y) ■■= D'^(x,y) « G^k^.^y) ° (4.4) 

where 



Likewise, the area element from S is pulled back to this on A by 

|X* A ^y\(^^y) dxdy = det (Gx|4,(^,j^)) det^ {D^i^^^y)) dxdy 
In particular, for any <I> E Diff{A), the surface area |S[ of S is given by 



|S| = 



j^det (^G'x|$(a;,j^)) det"^ {D^(^x^y)) dxdy . 



and is independent of Lemma l4. II below follows from definition: 



5 



Lemma 4.1. The diffeomorphism J^.^ is an isometry between S and Si if and only if 

D^ix,y) ° G^Ui.,y) ° = Gy|(,,j,) on A . (4.5) 

Likewise, li4-3{ ) is area preserving if and only if 

det (Gx|$(,,j,)) det" (D<i>(.,,)) = det (Gy|(,,,)) on A . (4.6) 

We now choose a particular reference parameterizations X, Y as follows. Recall that the 
Uniformisation Theorem [FK] implies that there exist conformal parameterizations of these 
surfaces. We may assume, therefore, that X (res. Y) are conformal parameterizations to S, 
(res. Si). This means 

G^ = fi(^l l) ' ^^ = "^(0 1) ^^-^^ 

where fi'^ (res. ry^) are the area densities on A associated with the conformal parameterizations 
of S (res. Si). With these particular parameterizations, condition (j4.5p for isometry of the 
mapping (j4.3|) is reduced to 

fimx,y)) mil^y^ = r?(x,y) J ° ) on A , (4.8) 

where |L'<I>p = o D*^. Likewise, (|4.3p is area preserving if and only if 

//($(x, y))det {D^i^^^y^) = r/(x, y) on A . 

Let us consider the set 

0(S, Si) :={<!>€ W/(A) ; ^($(x, y))det (i?<i>(,,,)) = r?(x, y) V(x,y)eA}. (4.9) 
Again, the result of Moser [M] implies that 0(S,Si) is not empty, provided 

|S| = |Si| . (4.10) 

Now, we recall (12. 4p . It implies 

/^(l?^(.,,)Gx|$(.,,)I)*<I>(.,,)) >0 (4.11) 
and, by Lemma |4. 11 it follows that S, Si are isometric if and only if 

k {D^i.,y) Gx!$(,.,,) ^*^(.,,)) ^ on A . (4.12) 
Under the assumption ()4.7p we obtain 

k[D<i>^,^y) G^\^(^,^y^D*<i>^,^y))=fi\^x,y))tr(^\D<^\l^^^^ . 
For $ e 0(S,Si) it follows 

GxU(,,,)I?*<I>(.,,)) =/^'($(x,y))tr -4ry2(x,y) . (4.13) 

Using (|4.11|4.12p and (f4T3]) we obtain 
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Lemma 4.2. Assume ^JTJW. //X,Y are conformal parameterizations and ^ € 0(S,Si) 
then 

H^{^) := I ^2($)tr (|L>^>p) dxdy > 4|Si| . (4.14) 
J A 

The equality in ^-14^ holds if and only if S and Si are isometric. 
We are now at a position to generalize (j2.8p : Calculate 



inf Ha('^) (4.15) 

and find the optimal mapping $ (if exists). 

As in section [3l the right action on the "manifold" 0(S, Si) is given by the Lie group of 
rj preserving diffeomorphisms on A: 

Or,{A) := {S G DiffiA) ■ rj{S{x,y))det (l)%j,)) = r?(a;,y) V (x,y) G A . } . (4.16) 

Lemma 4.3. The set 0^(A) is a group under compositions. For any $ G 0(S,Si) and 
S G 0^(A), <I> o 5 G 0(S,Si). In particular, Or,{A) acts on 0(S,Si) by right-translations. 

The Lie algebra associated with 0^(A) is 

Orj := {v G C7°°(A; M^) v • (r/^;) = on A , v h = on dA .} . (4.17) 

In complex notation ()3.4p , Lemma 13.11 implies that 

I? G if and only if v vi + iv2 = ir]~^dztp 

where ^p G C°°(A) and = on dA. 

As in the case of planar domains discussed in section [3l we consider now the flow G 
0(S,Si) generated by the Euler equation 

_ . = , G Or, , (4.18) 

and the Hopf function /$ defined, as in (jS.Sp . on A. The generalization of Theorem [1] is now 
formulated as 

Theorem 2. iei <I>° G 0(S,Si). Let 6e a flow i4.18\ ) where v'^^^ ~ irj^^djip^^^ satisfying 
^(t) ^ (^oo^^^)^ ^{t) =Q ondn for any t G M. T/ien 0(S, Si) for anyteR and 

li?^ (^$(*)) = 4am ^ ^2 f^^it)-^ (^-1^^) rf^dy . (4.19) 

Now, we generalize Definition 13.11 and Proposition 13. II as follows: 

Definition 4.1. A mapping Q : S ^ Si is quasi-rigid if and only if it can he decomposed as 

Q = Y*-^X-^ (4.20) 

where X : A ^ S, Y : A — > Si are conformal diffeomorphisms and <1> G 0^(A) verifying 
dHA{^)/dt = for any tp G C°°{A) for which tp = on dA. 

Proposition 4.1. Q given by ^.20^ is a quasi rigid deformation verifying ^f ^'^^ only 

if the corresponding $ G 0,y(A) verifies 

i)3md^{d^{fi\^)f^)r]~^} = on A. 



II 



i) 3m (/$§) =0 on dA. 
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5 Applications 

Theorem [1] in section [3] suggests an algorithm for calculating a quasi rigid deformations from 
a flat domain 17 C to another ili C where ()2.7p is assumed. 

Theorem 3. Let G 0(17). Define the flow ^>(*) in 0{n) given by the Euler equation 
hS. 'J\) for t > and i?^*) ~ id^tp^^^ G o(i7) is defined as follows: ip^^^ = tpo '^ + V'f^ + "^2*^ where 



V^(*) :=-am(9|/^(„) (5.1) 



on and 

as well as 
and 



#=0 , ^ = -:Jm(/,|) on 



then the flow 113. 3\) satisfies 



dt 



2 / / |jm(a|/<j,(t))|^(i2;dy 
Jn J 



an 



3m (4,) I 



2 



In particular, dH{^)/dt < and ^ is a steady state of ^3.3\) if and only if it is a quasi-rigid 
deformation. 

Similarly, Theorem [2] in section [J] suggests an algorithm for calculating a quasi rigid 
deformations from an embedded surface S C to another Si C M^: 

Theorem 4. Let S C M'^, Si C M'^ and a pair of confornial parameterizations X, Y verifying 

given by the Euler equation (J^J^ 
for t > and idzip^^^ G Ori is defined as in Theorem ??, where \5.1]} is replaced by 

# = -3mch{ch 7?-!} . (5.2) 

then the flow { 4-1^ satisfies 



dt 



ama,{a,(Ai^($W)/<,w)r?-i} 



dxdy- 



dA 



-V2($W) 



/n particular, dH{^)/dt < and $ is a steady state of (^TT^ i/ and on/y ifY^ ""^ 
quasi-rigid deformation ofT, into Si. 



2S a 



Remark 5.1. In case the surfaces S,Si are closed (e.g. both homeomorphic to the sphere 
S"^), then = ^q*'' as given in \5.^) . and there is no need of the component ipf^ , which 
take care of the boundary. 
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6 Proofs 

Proof, (of Theorem [1]) Consider an orbit induced by S^^^ G 0(17) := 0{Q; Q) via 

$W = $(0) o . 

A tangent of this orbit is given by the left representation 

S = voS ,ve o{n) , =1 . (6.1) 

Since V • = by definition it follows by (l22])that det{DS^*'^) = det (G^c*)) = 1 along this 
orbit. In particular. 

(y^<s>(o)os[']j ^ '^'st (G${o))(s(t)) • det {G s(t)) (^^-^ = 1 

so $W = o S-W e 0(Jl;Jli). In addition 

= [Dv oS\DS . 

In the left representation ^>(*+^) = o S'(^), 5W(x) = x + ri>(3;) + ©(r^), so 
We obtain from and ((631) 



(6.2) 



(6.3) 



+ T I tr 
In 



V dxdy + O(t^) . (6.4) 



However 



and integration by parts implies 



D 



D{tr 



vdxdy 



tr 



(V • v)dxdy = , (6.5) 



since v is divergence free. From (|6.4l 16. 5p we obtain 

-i^(^$W)= y tr (^Z)($(*))[L>^ + L>*i7]D*($W))(ixdy 



(6.6) 



Next, note that if Q is simply connected, then any smooth divergence free vector field can be 
written as 

V = v^ij := {-il^y,i^^), i^ec°°{n) ,ij = o on an. (6.7) 

It follows that 

^(Dv + D*v) = ^S{^j)+V{iP) (6.8) 



where, using □ := 3^ — 5^, 
S(V) := 

Note that, with 



nip 
□V- 



-'4'x,y 
V'x.s; 



SO 



while 



Hence, 



We now return to the complex notation. From (|3.41 13.5p . 

2 2 



xy 



hence, by 1^ and (|6l^[6l^ . 



Proof, (of Proposition I3.ip : 

We imply integration by parts on (j3.6p to obtain 



H ( = 4am 



d 

However, V = on dO, by 
d 



ipd^f^iDdxdy + 29^e / (V'^/${t) - Uwd^ip) dz 
in J / \Jdn 



so 



= 4Jm ^ y" il,dlf^(^)dxdy - 2^t j^^ f^it)Ch^dz . 



(6.9) 



(6.10) 
(6.11) 

(6.12) 
(6.13) 



□ 



(6.14) 



Part (i) follows from the first integral (j6.14p . since -0 is arbitrary in the interioir of i7. Since 
'0 = on 90, it follows by Corollary 13.11 that idzipdz is real valued on dfl. Hence 



(Udztpdz) = -3m (^h^idz^pdz^ = ±dntjj3m (^U^ 
Since dnip is arbitrary on dQ we obtain part (ii). 
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(6.15) 
□ 



Proof, (of Lemma 113]): Let S,Si G 0^(A). Then D{S o Si) = DS(^Si) ° DSi. Hence 
det {D{S o Si)) = det {DS)^^^^^ det (DSi). By definition, 

det{DS) = ri/r^{S) ; det {DS)(^Sr) = v{Si)/v{S o Si) 

so det {D{S o Si)) = ri/r]{S o ^i), so S o Si G 0^(A). The same argument holds for $ o 
where E 0(S, Si). □ 

Proof, (of Theorem [2]) : 

We repeat the proof of Theorem [3] up to (16. 4p . Here (16. 4p . is replaced by 
Ha (^$(*+^)) = Ha ($^*^) + ^ / (^^*^) (-D(^(*)) (I^t; + Z)*-?;) dxdy 

T {^l^l^ity) ■ v) tr (|l)cl>(*) dxdy + 0{t^) . (6.16) 
As in (j6.5p . the third term on the right of (j6.16p is reduced to 

j ;u2(^.W)D |tr ^ L'^>(*) I iJdxdy = - j tr L>^>W V • {n^{^^^^)v)dxdy . (6.17) 
Together with the forth term in (|6.16p we obtain 

(^$(*+^)) = Ha («>^*^) + ^ / (^^*^) (-D($(*)) (Z)^ + I?*t;) 1?*$^) dxdy 

- r ^ /^J$W)(V • v)tr (iDc&W [2) dxdy + 0{t^) 
Since € it follows from (j4.17p and (|6.18p 

(«>(*+^)) = -^A («'^*^) + ^ / (^^*^) (-C(^^*^) + Z^*^^) -D*$(*)) dxdy 

+ Ty" ^2($W)7?-i('j;-Vr/)tr (^|D$W[2)dxdy + 0(T2) 
We now use v ~ ir]^^dzip and (|6.7116.8p to obtain 



1?^; + D*v = r/-^S(V') + 2r/-^U('0) + 27?" 



Vxtpx-riy'4'y 
2 



(6.18) 



(6.19) 



By direct calculation 



2/i2 



„ „/, Vxll^x-Vyi^y 

'Ix Yy 2 



1^ 



W|2 



Vxi^X-Vy^y 

2 



y 
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while 



/x2($W)r?-^(t;.Vr?)tr(|Z)$W|2) =^2(^{t))^-2(^^^^_^^^^)(^|^ . (6.21) 

using (|6.10l6.1ip and (j6.18ll6.2T]) and differentiation at r = yield 



di 



dxdy 



(6.22) 

In complex notations (j3.4p . (j6.22p takes the form 



((DW) = 4Jm 1^ rj-^Udlijdxdy + j^t? (a^r?-i)(a^V)/*^ixdy 

= 4Jm y ^? (r/^^a^V) dxdy . (6.23) 



□ 



Proof, (of Proposition HTTj) : Integration by parts of ()4.19p yields: 



d_ 
di 



^a($(*^) =4am (^J^J ^W5^{^(^2(^(t))^_^^^^^ r^-^^dxdy 



2^e 



an 



V-'fi\'^^'^)Uwc^^^'Uz) . (6.24) 



The rest of the proof is, essentially, the same as that of Proposition 14.11 



□ 



Proof, (of Theorem [3] and Theorem Hj) : The proof of both Theorems follows from (16.14p and 
(j6.24p . respectively, where (j6.15p is applied in both cases. □ 
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